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Abstract. Employing the random matrix formulation of Chern-Simons the- 
ory on Seifert manifolds, we show how the Stieltjes-Wigert orthogonal poly- 
nomials are useful in exact computations in Chern-Simons matrix models. We 
construct a biorthogonal extension of the Stieltjes-Wigert polynomials, not 
available in the literature, necessary to study Chern-Simons matrix models 
when the geometry is a lens space. We also discuss several other results based 
on the properties of the polynomials: the equivalence between the Stieltjes- 
Wigert matrix model and the discrete model that appears in g-2D Yang-Mills 
and the relationship with Rogers-Szego polynomials and the corresponding 
equivalence with an unitary matrix model. Finally, we also give a detailed 
proof of a result that relates quantum dimensions with averages of Schur poly- 
nomials in the Stieltjes-Wigert ensemble. 



1. Introduction 

In the late eighties P, Witten considered a topological gauge theory for a con- 
nection on an arbitrary three-manifold Af, based on the Chern-Simons action: 

(1.1) Scs{A)^^[ Tr(AAdA+^AAAA A), 

47r Jm 3 

with k an integer number. One of the most important aspects of Chern-Simons 
theory is that it provides a physical approach to three dimensional topology. In 
particular, it gives three-manifold invariants and knot invariants. For example, the 
partition function, 



(1.2) Zk{M)^ J PAe^'^'^sCA)^ 

delivers a topological invariant of M, the so-called Reshetikhin-Turaev- Witten in- 
variant. Recent reviews are |21|3|. 

As reviewed in a great deal of interest has focused on the fact that Chern- 
Simons theory provides large N duals of topological strings. This connection be- 
tween Chern-Simons theory and topological strings was already pointed out by 
Witten |S] (see also 0), and then extended in 0. 

Recent progress in Chern-Simons theory includes a description of Chern-Simons 
theory on certain geometries in terms of models of random matrices. Consider the 
partition function of Chern-Simons theory on a Seifert space M — X{^, . . . j^-)- 
This is obtained by doing surgery on a link in with n + 1 components, out of 
which n are parallel, unlinked unknots, and one has link number 1 with each of the 
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n unknots. The surgery data are Pj /qj for the unhnked unknots, j = 1, . . . , n, and 
for the last component. The partitfon function ifQ 4j: 

(1 3) Zcs(M) ^ ( ^] NS"(^)]'"^'e^.isn(./P)-^ 



121 



E 



This expression gives the contribution of the reducible flat connections to the par- 
tition functions. Recall that for both 5''^ and lens spaces this amounts to the exact 
partition function. The case n — corresponds to the three-sphere that leads 
to (|1.5|l . Thus, for the case of U{N), and focusing on a particular sector of flat 
connections, we get the following matrix model: 



(1.4) 



N OO n" n 2sinh 

iV- n.<.(2sinh^:" 



Of course, the simplest case is that of with gauge group U{N), which is given 
by the partition function of the following random matrix model: 

(1.5) z=L^i^/ne-"'-n 

1=1 i<j ^ ' 

From the point of view of topological strings, this describes open topological A 
strings on T*5^ with N branes wrapping 0]. This latter shown in [5], 

can be studied with usual techniques of random matrix theory. More precisely, the 
Stieltjes-Wigert polynomials, a member of the g-deformed orthogonal polynomials 
family jU], allows to compute, in exact fashion, quantities associated to the matrix 
model. In the computation, the g-parameter of the polynomials turns out to be 
naturally identified with the g-parameter of the quantum group invariants associ- 
ated to the Chern-Simons theory. This is so because the previous model can be 
easily mapped into: 



(1.6) Z = y'[dM]e-^T^('°s^^)\ 

named Stieltjes-Wigert ensembles, after the associated orthogonal polynomials. 

Chern-Simons matrix models have been further considered in {10^ and ll] -|17 | 
and also play a central role in q-2V) Yang-Mills theory 5H|-|2SI- Most of these 
works focus on the relevance to topological strings. In jHIEIj, the emphasis is on 
exact solutions and on the special features of the matrix models. The works of 
Caporaso et al. [201 also make an extensive use of the properties of the 

Stieltjes-Wigert orthogonal polynomials. We shall be focussing here on aspects of 
the Chern-Simons matrix models that have to do with the associated system of 
orthogonal polynomials. 

This paper is organized as follows. In the next section we shall construct a 
biorthogonal extension of the Stieltjes-Wigert polynomials, in order to study the 
matrix model (|1.4|l when n — 1 and n — 2. These polynomials have not been 
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discussed in the (vast) orthogonal polynomials literature, so most of our effort is 
on their derivation and to establish some of its fundamental properties. They are 
necessary if one wants to obtain full analytic results when the geometry is some- 
thing more complicated than S^. Note that matrix models in the lens space case 
have already been studied (with loop equations) |13| , but if one desires an all order 
result as in 8 , the knowledge of orthogonal polynomials is then necessary. After 
the construction of the biorthogonal Stieltjes-Wigert polynomials in Section 2, we 
discuss some of their mathematical properties in Section 3. In the last Section, 
we discuss several aspects of the Chern-Simons matrix models by focussing exclu- 
sively on properties of the (ordinary) Stieltjes-Wigert polynomials. In particular, 
we clearly establish the relationship with the discrete matrix model that also ap- 
pears in g-2D Yang-Mills theory, and also employ the intimate relationship between 
Stieltjes-Wigert and Rogers-Szego polynomials to find the exact relation between 
the Stieltjes-Wigert matrix model and Okuda's unitary matrix model ^H]. Finally, 
we give a detailed proof, employing a mixture of combinatorial and orthogonal poly- 
nomials results, of the equality between quantum dimensions and averages of Schur 
polynomials in the Stieltjes-Wigert ensemble 0. We conclude with a summary and 
with some avenues for further research, presented in the Conclusions and Outlook. 



2. Biorthogonal Stieltjes-Wigert 

Let us consider the generic expression (jl.4(l in the n — 1 and n — 2 cases, that 
correspond to the case of lens spaces. We are lead to a biorthogonal extension of 
the model: 

(2.1) ../n.-..f/-n(^^i"..^)(-'.Hi^^)$- 

Recall that a biorthogonal ensemble of random matrices has the probability density 
El: 



N 



(2.2) P(xi, xn) ^CnI[u; (x,) [] {x, - x,) (xf - 4) , 



where fc is a fixed real number. In total analogy with the usual Hermitian case 
(fc = 1 ) one can study H2.2|l by considering a pair of biorthogonal polynomials: 



(2.3) J LO (x) Yn {x, k) Zm {X, k) dx = hn,,kSn,' 

with: 

(2.4) I Yn{x,k)x''^L0{x)dx = a^n^5n,j. 

Zn {x, k) X-' U! (x) dx — Pl^^Snj- 



We warn the reader that the term biorthogonal is employed in different contexts 
in the literature. The classical cases (Hermite, Laguerre and Jacobi) were worked 
out in Note that (|2.2|l is exactly the type of ensemble that ()2.1(l leads us to 
consider since: 
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i i<j ^ 

(2.5) = .-"^/n|^^-^^"^^^'n(^/.^'"-2/r)(2/..^'''-^r)' 

i i<j 

with Ui = log Hi, = l/2gs and a = -1 - , /3 = ;^ + ^- Finally, with 

= e^K^p and some rewriting: 

(2.6) 

which is of the form H2.2|l with the log-normal (Stieltjes-Wigert) weight function: 
Lo {x) = e"" ^ the g-parameter is then q — e~ i^-^p'-^ = . Therefore, if we 

want to go beyond the 5*^ case one has to construct the biorthogonal Stieltjes-Wigert 
polynomials, not available in the literature. Thus, this is our main task in what 
follows. The method we have chosen is based on a simple but fundamental result by 
Askey, that relates the g-Laguerre orthogonal polynomials and the Stieltjes-Wigert 
polynomials |25| : 



(2.7) \my Ll{q-^x-q)=Sn(x-q), 

and then we take into account the biorthogonal construction of the o-Laguerre poly- 
nomials, carried out by Al-Salam and Verma in the early eighties Q The Stieltjes- 
Wigert polynomials are 27^ 

n 

(2.8) S^{x\q)^——Y, 

The limit (|2.7|) will provide us with a biorthogonal extension of the SW polyno- 
mials starting with the g-Konhauser polynomials. Therefore, following [20], let us 
write: 



(2.9) {x, k\q) EE \ y ^-^ '\>i'^ ,^ r 1^ 1 ^ ' > 

and 



(2.10) ri")(x, k\q) ^ytT. ^r, 

Wn ^ [q]r 

with 



The resulting polynomials were named g-Konhauser as they could also be interpreted as a 
(j-deformed version of the biorthogonal Laguerre polynomials, worked out by Konhauser. 
^In 1251 they appear, in Eq. (2.5), slightly reformulated. 
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(2.11) 



s=0 



These polynomials satisfy: 
(2.12) 



with fcl") 



1+al „ — nk 



[<l]n 



with respect to the normalized q-Laguerre measure. We have to study: 

(2.13) Zn{x,k\<l)= lim 4")(g-"^,A:|g), 

Yr,{xM<l)= lim Yi^\q-"x,k\q). 

a — ^oo 

In the first case, one readily finds: 



1 " (^q-nk. qkyq^kj(kj-l)+k](n+l) 



(2.14) 

which can be conveniently reexpressed 



,kj 



(2.15) 



1 " 



Regarding Yn{x, k\q) we have to find 6,- = lima^oo q "'^fc". Employing Taylor [221 
one can write: 



(2.16) 
therefore: 



.=0 t'?]^ 



(2.17) 



Taking the a — > oo limit and using the finite g-binomial theorem [S]: 



(2.18) 
so that 



x"^ 

[qx; q% = ^ - ^^(-l) 



^|fcr(r-l)+r^r 



(2.19) 



(-1)^ 



■|fcr(r — l)+r 



For later use note that one also have: 



6 



YACINE DOLIVET AND MIGUEL TIERZ 



(2.20) 

From this one gets: 



1 / d 



(r) 



(2.21) 



Y,,{x,k\q) = —J2{~1Y 



^^r(r+l) + ifcr(r-l)^r^ 



For k = 1, both polynomials reduce to the Stieltjes-Wigert polynomials H2.8(l . 
Writing Yn{x, k\q) = yn,kX^ + ■■■ and Zn{x, k\q) = Zn,k x^'^ + ... one finds: 



(2-22) Zn,k = q 



( -nk. k\ ( -i\n in^fc(fc+l) 

^ 19 ,g J« ^lfcK(fc«-l) + fcrt(n+l) ^ g3 



and 

(2.23) ^ 
This leads to: 



1 



n{n-l) " \ I I -'-J g'' 

n\[q\n \dx J [q\n 



(2.24) 



< Yn{x,k\q),Zm{x,k\q) >= /i„(5„,™, 

Ada; 



with respect to the measure t — , , , , — , with A such that < 1, 1 >= 1 and 

g — nk 



(2.25) 



[q]n 



Using this, we can find for example: 



(2.26) Z^^'? = 7V!(^)^/V^'"^'+^^'+^^^'^"'^^'+'+^^'^'~'^^ 17(1-'?*)'^"^ 
27r -'- -'- 

which reduces to the known formula when P = Q = 1 011]. 

3. Mathematical properties of the biorthogonal polynomials 

Since the biorthogonal Stieltjes-Wigert polynomials have not been addressed in 
the literature, we derive here some of its fundamental properties. 

3.1. Behavior under dilatation. First, we find some generating functions for 
Zn{x; k\q) and y„(x; k\q) (t ^ q-^) 



(3.1) 

and 



Y^Znix-Mqy 



ri>0 



fjtx'^) 

{t;q'')oo 



(3.2) 
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with 
(3.3) 



r>0 '3 r>0 



dx 



(r) 



{qx\q^)n\x=o- 



We rely for this on formula (4.2) from [2^1 ■ The expression for Z is essen- 
tially property (4.1) in j^. For (|3.2I) we use the explicit expression obtained for 
Yn{x]k\q). Let us introduce the moment generating function: 

(3.4) G(i,x)EE^ J^r„(x;fc|g)r, 

that can be written as: 



(3.5) G{t,x)^Yl 

r>0 



q 2 




r! 




■(r-l) 


q 2 




r! 




-(r-l) 


q 2 



n'>r 



{q^]q'')r 



x=0 



r>0 



dx 

2 



E 

n>0 



('■) 



(i;<z'=)o 



-(xt)' 



rfa; 



(r) 



(ga;;q'=)„|:r=o ; t^q 



-k 



In the second line, the extra piece we add, being a degree r — 1 polynomial in x 
does not contribute due to the derivative. In the third line we use the g-binomial 
theorem, and in the fourth one we make the change of variable x — » xt. 

Now, by taking 1)3. 2|l with x Xx, introducing in the r.h.s the factor 
and matching the coefficients of on both sides one gets: 



(t;9'=)o 



(3.6) 
with: 



y„(Ax;fc|g) = ^7„^.(A)y,(x;fc|g), 



(3.7) 



7„,-(A) 



[q], (gfc;gfc)„ A-'(A;g^)„_, 
[q]n {q^]q^)j {q^;q'')n-j 



and similar steps involving Ea. l3.1l give: 



(3.8) 
with 



Zn{Xx]k\q) = ^Q^j{\)Zj{x;k\q), 



8 



YACINE DOLIVET AND MIGUEL TIERZ 



(3.9) 



C„,(A) 



One has C- 



71J 



7„j for fc = 1 as it should. Moreover 7„ 



A", by 



matching the dominant coefhcients on both sides. Note that equation (4.2) in [2^] 
contains a typo as it does not fulfill this last condition (it would give C„„ = !)• 
And of course one has Cnj(l) = 7ry (l) = ^nj- 

Even though and are defined for j < n we extend for convenience their 
definition through 



(3.10) =7„^. =0 if j>n. 

3.2. Recurrence formulae. The Stieltjes-Wigert polynomials (|2.8(l satisfy: 



(3.11) 



Sn-iix\q) - (1 - q")5„(x|q) + xq^' Sn-i{xq\q), 



an identity used by Chihara in [28j . to prove that the zeros of the polynomials 
satisfy: 



(3.12) 



*^n,m ^ Xn — l^m ^ qXn,m+l^ 



where n denotes the order of the polynomial and m indexes the zero. Note that the 
zeros of the SW polynomials are an interesting quantity in the context of topolog- 
ical strings In what follows, we find the same identities for the biorthogonal 
polynomials. 

3.2.1. Fundamental recurrence relation. Note that for the particular value A = q^^ 
(123 gives: 

(3.13) C„,('Z"')-0 if j<n-2. 

This implies the following simple recurrence relation for Zn{x, k\q): 



(3.14) 



Zn{x, k\q) - Zn-i{x, k\q) = q^'^Zniq ^x, k\q) 



Certainly, if one writes the Zn{x,k\q) — '^j^qTtijx''^ , it can be checked directly, 
from the explicit expression in 12.15|1 , that one has: 



q ' r, 



(3.15) ''"nj' ''"n— l,j 

which implies the recurrence relatior^. 

For fc = 1, H3.14|l reduces to the following relation for the Stieltjes-Wigert poly- 
nomials 



(3.16) 



5„(2/)-5„_i(y)=<z"5„(g-iy). 



''incidentally, this a check that Eq. 13.91 is correct. 
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3.2.2. Moment generating recurrences. From H2.19|l one has {bn,o — 1, &n,-i = 0): 



(3.17) 



^n+l,r ^n.r nk-\-l ^n,r — l 



[q]r [q]r [q]r^i ' 

which imphes the following recurrence relation for the Yn{x,k\q) : 



(3.18) (1 - g"+i)r„+i(a;, k\q) = Y,,{x, k\q) - g"'=+ixy„(gx, k\q), 

or, equivalently: 



(3.19) 



xYnix,k\q) - q-"'' {Y^{q-'x,k\q) - (1 - q''+')Y^+,{q-' x, k\q)) 



We proceed in analogous way for Zn{x; k\q). For convenience we introduce coef- 
ficients Cn,r, such thal|j: 



(3.20) 
that is: 



Z4x,k\q)= 



^ ^" ^ j,kr q^kr(ki — 1) 



(3.21) 



Then, as in the previous case c„.o 1, c„^_i = 0, then: 



22) C.n+l,r _ r ^ ^nk+k '^n,r-l 

{q'';q'')r {q'';q'')r {q'';q'')r-i'' 
and one gets the recurrence relation for Zn{x; k\q) : 



(3.23) (1 - g'^("+i))Z„+i(a;,A:|g) - Zn{x,k\q) - g"'=+^a;'=Z„(g'=x, A:|g), 
equivalently: 



(3 24) 

x^Z^ix, k\q) = [Zniq-^x, k\q) - (1 - g'^'^+l' )Z„+i (g-'^x, k\q)) . 

One can easily check that these recurrence relations both reduce, taking fc = 1, to 

To conclude this Section, since we know the explicit behavior of the polynomials 
under dilatation, we can employ (|3.7|l and (|3.9(l . and then, H3.19|l and (|3.24|l , to 
obtain an explicit way to compute the moments < x''Yn{x, k\q)Zm{x, k\q) >. For 
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instance, one has: 
(3.25) 

< xY^{x, k\q)Z^{x, k\q) > - q-^^'^k^ {in.^iq-') - (1 - 9"+')7„+i,™(9~')) 

„— (n+m)fc 

= hn.rniq-') " (1 " 9"+' )7„+l,™ (<Z"' ) • 

Note that this does not work so well when A; = 1 as 7„ mi'l^^) well defined 

in this case, according to the discussion above. 



4. Stieltjes-Wigert, other results 

4.1. Moment problem and q-2T) Yang-Mills. The Stieltjes-Wigert matrix model 
posses distinctive mathematical features, in comparison with other, more usual 
models in the literature, such as matrix models with Gaussian or polynomial po- 
tentials. The log-normal weight function leads to an indeterminate moment prob- 
lem |3(JI nT\ and consequently, the Stieltjes-Wigert polynomials are not dense in 
L'^{x,w{x)) (see [HI El for details). One of the consequences, discussed in 10 , is 
the discretization of the Chern-Simons matrix model. This opens the possibility for 
the discrete and continuous versions of the same model to share the same orthogo- 
nal polynomials. Note that in '21', we find the suggestion of studying the discrete 
matrix model (see below) with orthogonal polynomials, as done in (32j for the 2D 
Yang-Mills theory case, that lead to a discrete Gaussian matrix model. Recall that 
Gross and Matytsin '32' found a discrete version of the ordinary Gaussian matrix 
model in their study of the 1 /N expansion of the partition function of 2D QCD on 
the sphere: 



-t-oo 

(4.1) Z{A,N)= e-^^^^"Y[{u,-Ukf. 

ni,...,U7v — — oo j<k 

In contrast to the continuum case, that is solved with Hermitc polynomials, the 
discrete Gaussian weight does not have a closed system of orthogonal polynomials 
associated (see [33], for a recent discussion of discrete matrix models). So, they 
could not rely on known orthogonal polynomials, hence the difficulty of studying 
the discrete Gaussian model. Actually, the large N phase transition of the theory 
is related with the discrepancy between the discrete and continuous orthogonal 
polynomials. 

The orthogonal polynomials for the discrete matrix model in the Chern-Simons 
case are the Stieltjes-Wigert polynomials, as we show below. Note that the model 
in Chern-Simons/g-2D Yang-Mills differs from the Gaussian/2D Yang-Mills model 
in the exponentiation of the eigenvalue repulsion at large distances, {uj — Ufe)^ — > 
sinh^ (ui — Uj) . But this exponentiation is precisely the ultimate responsible of this 
continuum/discrete equivalence (see [TU]). 

The following detailed computation highlights this special property of the Chern- 
Simons matrix model. Let us proceed then to show the details of the derivation 
from Eq. (25) to Eq. (26) in ^ 
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+00 

(4.2) Za= e-*S.«? J|4sinh2(|(u,-z.fc)^ 

+ 00 

ui,...,un — — oo i<fc 
+00 

«!,...,«« = — 00 j<fe 

+00 N 

= E n 9^"'+"' n('?"^' 

ui i.-.j-fiAf — — 00 z— 1 j<k 
N +00 



\2 



i— 1 n— — cxD j<k 



N 



•^0 .=1 i<fc 

where we have introduced c = e^^° = and we recah that we have as usual 
q = e~^= . Moreover, one has: 



+00 



(4.3) ^,(^)^_1_ J2 c-q'^+^Six-cq"), 



n— — 00 



which is the discrete measure equivalent to the continuous distribution w{x) in 

as far as the integer moments are concerned |28l 1341 . The normalization is given 

by: 

(4.4) M(c) ^ (-cg3/2, _c-lg-l/2, g; q)^ = hcq'^'U[-c-\-'^'U[qU 
therefore, we can use the equivalence between these two measures to write: 

(4.5) Zd = c^(i-^)gf Af(c)^ / l[dx^ w{x.O [] (x, - Xkf . 
However recall that, in the simplest case P = Q = 1 in (|2.6(l , one has: 

(4.6) z''' =q-'^ i^^y l[dx,w{x,) Hix.^Xkf. 

•'^ 1=1 j<k 

Thus, one has a quite simple relation between the discrete and continuous Stieltjes- 
Wigert ensembles: 



^Note that even though there are infinitely many discrete measures ui^ equivalent to ui (the 
parameter c is a real number), in the Chern-Simons case we are discussing, the constant c is a 
function of . 
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(4.7) 



-— e 29. 



27r 



I 2sinh(- 

-i(l-2Af+37V^) 



E 



Note the inversion of the couphng constant between the l.h.s. and r.h.s 



e~*i:.»? J](2sinh(|(n, -n,.) 

ni,...,njv = — oo j<fe 



4.2. Prom Stieltjes-Wigert to Rogers-Szego: unitary matrix model. The 

Stieltjes-Wigert polynomials turn out to be intimately related to the Rogers-Szego 
polynomials I^B! , that are orthogonal on the unit circle. This is useful to estab- 
lish in detail the exact relationship between the Stieltjes-Wigert matrix model and 
the unitary model considered by Okuda 15 . Following |37| . we recall the defini- 
tion and relations between the Rogers-Szego and Stieltjes-Wigert polynomial. The 
Rogers-Szego polynomials are defined as: 



(4.8) HMq)^Yl I 

and they satisfy an orthogonality relation on the complex unit circle: 
(4.9) 



J\w\=i \ 2i 



w 



where 03{z\q) is the third Jacobi theta function (see definitions in the Appen- 
dix). Note that the orthogonality coefficients hm = -^^k^ are identical to the ones 
(Stieltjes-Wigert) that directly give the Chern-Simons partition function in the S'^ 
U{N) case This is enough to write down an unitary matrix model for the 
Chern-Simons partition function. However, let us show this point with detail. The 
polynomials are also orthogonal with respect to a measure defined on the full real 
line EH: 



Mr 



1 /"TOO 

(4.f0) / i7„,(-g-i/2e-2»^- | q)H,,{~q-'^'e'^^- \ q)e-- dx 
introducing through q = . Now consider the Stieltjes-Wigert polynomials 

Sn{x) 27 J 



(-i)"ff»/2+3 " r 
(4.11) S^ix)= ^ ^5- 

viiU 



(-l)"g"/2+3 
\/\<lU 



Sni-\/qx\q), with Sniz\q) = ^ 



fc=0 



These polynomials fulfill the following orthogonality relation on the real line: 



CHERN-SIMONS THEORY AND STIELTJES-WIGERT POLYNOMIALS 



13 



(4.12) ^ S^i-q-'/'e-'^^ I q)S^i-q-'/'e~'^^ I q)e'^'dx = 



Using an elementary property of the g-binomial coefficients, tlie two equivalent 
relationship follow: 

(4.13) Hn(x\q-^) ^ Sniq-^'xlq) and ^„(x|g-i) = i/„(q-"a;|g). 
Then: 

(4.14) 

k 



< Pn,Prn > u, = / e"^ (-gi/^z)^,, (-gi/2^)dz 







/ e-'=^(--^)^S„(-<Z^/^c-)^n(-g^/V)da 



1 n-\-OQ 

p^^^q-^l^— \ e--'^„(-g-i/2e^)^„(-g-i/2g^)da: 

VT^ J-oo 

/ e-^^i?„(-g-V2e2'A<-)i/„(_q-V2e-2.M.)dy 

^ ^_^H™(-g-^/^T^; I g)iJ„(-g-i/^^ | 9)63 (^IV? 



where we also have used, between lines 3 and 4, the fact that Sn{ae^^'^^ \ q) and 
Hniae'^'^'^y I q) are related by a Fourier transform We also introduced 2/i = i 
and: 

(4.15) P™.„ = (-1)™+"^==^. 

The next line is given by the results of the previous section. Now consider 
Eq. (3.22) in [ini, it reads: 

(4.16) 

N 



\i=l / i<j i<j 

I (feQoo(e-1,)) ^^det^(il,_,(e-0)^^det,(i^.-i(e--)), 



\W\ 

N(N-l) 



where 



(4.17) eoo{e^'\q) = J2l' 



e ■' . 
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Then, considering that Jl] and PZ] have different conventions for the third Jacobi 
function one sees that: 

(4.18) e(?)(e'«k) = e(^)(^IV^). 

Therefore, one can continue the computation and write: 

(4.19) zcs=^-^^ <[^jf^^{{~iy-\-^^,p,^M))) >», 

which then connects with the usual expression of the partition function in terms of 
the orthogonal polynomials for the measure on the real line. 

Incidentally, both Stieltjes-Wigert and Rogers-Szego can be interpreted as the 
ground-state wavefunction of a (/-deformed harmonic oscillator '37' . This is an ap- 
pealing property as it has been recently shown that the Stieltjes-Wigert polynomial 
describes _B-brane amplitudes on the conifold |29|. 



4.3. Quantum dimensions as averages of Schur polynomials in the Stieltjes- 
Wigert ensemble. In this section we prove a formula for the averages of Schur 
polynomials that appears in pil39|. without relying on the equivalence with Chern- 
Simons theory. 

Recall that Schur polynomials Sa ,38 constitute a basis of symmetric functions 
in a given set of variables x = {xi) and are indexed by Young diagrams A. If the 
variables x are seen as eigenvalues of some matrix M £ sin then 5\{M) = TrA(M) 
is the trace of M in the representation associated to A. The Schur polynomials 
may also be more directly defined in terms of the skew-symmetric polynomials 

= det(x^^ " "*) as: 

(4.20) .a(.) ^ 

where as{x) is the Vandermonde in the variables x. The result we want to show is 
the following: 



(4.21) < Sa(A/) g-«l^l-^<'"'2?;,. 
with 

(4.22) C,^("' = {n + 1)|A| + ^^(A? - 2^K), 

i 

the Casimir of the representation labeled by the Young diagram A and |A| its total 
number of boxes. Background material for this Section is presented in Appendix 
IbI Note that this quantity can be rewritten using Ea. l|B.6ll as: 

(4.23) C^^"^ = n|A| + 2(n(A') - n(A)), 
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where A' denotes the conjugate partition. The quantum dimension is defined by 
the g-hook formul^l 



(4.24) 



n 



\n + c{x)\ 
[h{x)\ 



where for each box x — {i,j) of the diagram h{x) = Xi + X'j — i — j + 1 is the 
hook-length and c{x) = j ~ i the content of x. 

4.3.1. Case of 1- column diagrams. The quantum dimension of the j-th fundamental 
representation of sZ„, which is associatecQ to the partition (1-'), or a one-column 
Young tableau of length j, is 



(4.25) =dim,A(,) = 

Moreover, the monic Stieltjes-Wigert polynomials can be writter0: 



(4.26) TTn{x) = ^(-l)"-J(;0-")(j+"+i; 



x^ =< det(x - M) >i 



Besides, the following formula holds for the characteristic polynomial: 



det(a; - M) = ^(-l)"-^5(i„-,)(M)x^ 



(4.27) 

with 5\{M) the Schur polynomial associated to the partition A. Therefore: 



(4.28) Y.^-!)'"-' < S(i.-.)(Af) >^ = ^(-l)"-JgO"-«)0+«+5) 



J=0 

from which we extract: 



X-' , 



(4.29) 


<5(p)(M) >„=q-^(2„-, + i) 


n 




n 






J. 


q 


J. 


1 



Using l|4.22|l . one sees that Eq. H4.29(l is indeed consistent with Eq. I|4.21(l . 

4.3.2. General case. To study the case of a general Young diagram we note that 
as a generalization of Eq. 14.271 higher powers of the characteristic polynomial are 
generating functions for diagrams with a higher number of columns. Relying on a 
formula first computed in 41 we then relate the average of Schur polynomials to 
some determinant of Stieltjes-Wigert polynomials. From (1.4 Example 5 p. 67) 
we see (taking a slightly more convenient notation) 



^See for instance §4.4 in (4Q| for a clear discussion of the definition and properties of quantum 
dimensions. 

®We will freely switch notations between Young diagrams and partitions in the following. 
^The measure being here normalized such that < 1 >w= 1- 
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(4.30) 



i=lj = l A;Ai<fc 



where A = (Ai, . . . , A„) is a Young diagram with at most k columns as imposed by 

the condition Ai < k. The associate diagram A is defined as (n — A„, . . . , n — Ai). 
Therefore, if one considers —yj to be the eigenvalues of M, one immediately gets: 



(4.31) 



[|det(x,-Af) - J2 (-l)'^'sA(M)s3,'(a 

i=l \-\i<k 



By a standard result on characteristic polynomials j41| . we have: 

7r„(xi) ... Tr„+k-i{xi) 



(4.32) < TTdet(x, - M) 



7r„(xfe) ... TTn+k-lixk) 



with as{x) the Vandermonde determinant of the x variables. 

We now turn to the r.h.s., that we call A for convenience. From the explicit 
expression for the Stieltjes-Wigert polynomials one obtains: 



-ij(2«+2<T0)-2-ij + i 



ii,...,ik cr&&k 



n + cr(j) - 1" 



ri+o-(j)-l-ij 



io-e6fc j=i 



ij + o-(j) - 1. 



n 

i=i 



Zl,...,4fc 



X det 

l<a,h<fe 



n + b-l 
ia + b-l 



A; 

i=i 



in the second line we have relabeled ij ij + cr{j) ~ 1. Now we study the deter- 
minant in the previous expression and show that if ii > . . . > i„ : 



(4.33) 





'n + b-1' 


) = A (A) 


n 




det ( 

l<a,b<k \ 


ia + b- 1 


A 


1 
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with the constant Aq{X) = g"^'^ \ as we show in appendix IbI Here A' is the partition 
conjugate to A and equal to (ii, 12 + 1, . . . , «fc + fc — 1). [^] is a notation generalizing 
the g-binomial coefficients ["] and which is defined by the g-hook formulJ^ 



(4.34) 



1 _ qn+c{x) 

xex 



This is nothing but the analog of quantum dimension Eq. (|4.24f) where instead of 
using the [.J version of the g-integers one rather use0 [.]. Identifying with the 
l.h.s. of Eq. ()4.31|l we obtain: 



(4.35) 



The last step we need to perform now is to convert [^] in terms of T^x and check 
that the prefactor is given by Eq. 14.221) . To this end we ffi-st note that due to 
Eq. l|R8l) and Eq. l|B7jl we have: 



i(n~l)[A|-n(A) 



(4.36) , = q 
and 

(4.37) ^ ^i,{2n-i, + i) - (j - l)(2n + j - i) = -(2n - ^)|A| + ^ A;^ - 2jA;.. 

j 3 
To rewrite things in terms of the partition itself, rather than its transposed, we use 
the relationship: 

(4.38) Xl - 2i\, = 2{n{X') - n(A)) - |A|. 

i 

Collecting all the prefactors we can eventually write our final result: 
(4.39) 



which coincides with Eq. (|4.21|) . 

5. Conclusions and Outlook 

We have constructed the biorthogonal Stieltjes-Wigert polynomials, necessary 
for computing expressions such as H2.1|l . that appear in Chern-Simons theory. The 
polynomials are not discussed in the mathematics literature, so a great deal of the 
effort has been devoted to the explicit description of their fundamental properties. 

The construction of the biorthogonal Stieltjes-Wigert polynomials may very well 
be a necessary technical step for the computation of knot invariants in exact fash- 
ion employing orthogonal polynomials. Note that, so far, the topological invariants 
computed with orthogonal polynomials are only Chern-Simons partition functions 
(more precisely, only the case of with gauge group U{N) "E). Indeed, accord- 
ing to Marifio 021 1 the results in 0j can be extended in order to obtain random 

"'^'^We warn the reader that for convenience we adopt a slightly different notation for [^] com- 
pared to 1381 in the sense that its value for the partition (f) is the usual Gaussian polynomial 
whereas in EHl, = [-^]. 

"'^"'^ Recall that [n] = q ^ [nj 
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matrix descriptions of other Chern-Simons observables. The case of torus knots for 
example, amounts to: 



-(pQ(A^-p^) + (f + ^)p^ 



\PQ\^m 



i—1 i<j ^ ^ ^ 

where S\{xi) are Schur polynomials associated to the partition A (representations 
of U{N) are labelled by partitions A). 

That is to say, the 1)2. l|l solved here, but with an insertion of a Schur polyno- 
mial. However, an obstacle can be the lack of a computational device for random 
matrix-like quantities with such a term. Note that the case of an ordinary Gaussian 
Hermitian matrix model with a Schur polynomial insertion, was solved in |43| by 
purely combinatorial methods, with no use of Hermite polynomials at all. Never- 
theless, as we have seen in the last Section, a computation of the Stieltjes-Wigert 
ensemble with the Schur polynomial can be carried out with a mixture of combi- 
natorial and orthogonal polynomials techniques. Therefore, it turns out that we 
have studied in detail the two cases comprised in (|5.1() . The biorthogonal case with- 
out the Schur insertion and the average of the Schur polynomial in the orthogonal 
(P = Q = 1) ensemble. 

We have also studied other aspects of the (ordinary) Stieltjes-Wigert polynomials 
that are of direct relevance for the corresponding matrix model. In particular, 
we have discussed in detail the equivalence of the Chern-Simons matrix model 
with its discrete counterpart, of very much interest in q-2D Yang-Mills theory, 
and also discussed the close ties with Rogers-Szego polynomials, which are defined 
on the unit circle (both sets of polynomials being an equivalent solution of a q- 
deformed harmonic oscillator problem) . This relationship allows to clearly establish 
the relationship with the unitary matrix model discussed in (T5) . Fundamental 
properties of the Stieltjes-Wigert polynomials like their asymptotic behavior and 
the above mentioned q-deformed harmonic oscillator property, may be of interest 
in connection with the recently established role of the polynomial in the study of 
topological strings . We hope to address some of these issues in future work. 

Acknov^rledgments 

We arc grateful to Edouard Brczin, Bernard Julia, Vladimir Kazakov and Marcos 
Marino for valuable comments and correspondence. 



CHERN-SIMONS THEORY AND STIELTJES-WIGERT POLYNOMIALS 



19 



Appendix A. Normalizations 

To follow standard conventions it is convenient to have the orthogonal polyno- 
mials either monic or normalized to unity. Hence, we rewrite the previous ones a 
little bit (it will also make the link with the usual Stieltjes-Wigert for fc = 1 more 
transparent). 

A.l. Notations, k = 1. From Szego we have for the Stieltjes-Wigert polyno- 
mials 



(A.l) 



(^l)"g"/2+l ^ 
Pnix) = 7= 2^ 



with ["] the g-binomial coefficient: 



[<l]r 



These polynomials are orthonormal for the scalar product <, >u, induced by: 



(A.2) ^^a:) ^ ^ e-^ \ 

with q = e~^/^'' as usual. Note that one has: 



(A.3) <1,1>»=1/V9- 

The Sn polynomials in Eq. (|2.8() to which Askey refers as the SW polynomials 
are written in a slightly different form. They satisfy 

(A. 4) < Sji: Sm "T^j ^n,mj 

with < , > the scalar product associated to the measure "^['j^g/^; ] — with A a 
normalization constant such that < 1,1 >= 1. 
Then, the polynomials defined by: 

(A.5) Sn{x) = (-l)"v^g"/'^n(x), 

are orthonormal for <, >. One then sees that: 

(A.6) Pn{x)^q^/^Sn{Vqx). 

Since Al-Salam and Verma "^Sj have the same notations as Askey for the g-Laguerre 
polynomials, we will make the same rewriting when using the biorthogonal Stieltjes- 
Wigert polynomials in the context of Chern-Simons theory computations. 

A.2. k arbitrary. Normalizing and changing variables as in the previous section 
we define new polynomials 



(A.7) 
and 



(_l)"gl/4 

Rn{x, k\q) = — — Yn{^x, k\q) = rn,kx" 
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(A.8) 

Then one has 



T„{x, k\q) = ^ 'J. ZniVqx, k\q) - t^^kx''" 



(A.9) 

and 



rn,k 



(A.IO) 



which reduce to - — , when fc = 1 as expected. Then one has 



(A.ll) 



< Rn{x,k\q)Tn{x,k\q) 



FoUowing Eq. H3.19|) and H3.24|l the recurrence relations for these orthonormal poly- 
nomials read 

(A.12) 



xRnix,k\q)=q-"''-'/^ (^Rniq-'x, k\q) + g^'/Vl " q"+'Rn+liq-'x, k\q)^ , 



(A.13) 



x''T^{x,k\q) = lT,,{q-'^x, k\q) + q-'^/' {1 - ) y^T„+i (g-'^a;, fc|g) 



Appendix B. Proof of Eq. (|4.33|) 
Eq. (|4.33ll is not obvious at first sight because when q = I, Giambelli's formulJ^ 



would lead us to write: 



(B.l) 

However, classically one also has 



det , 11= diniA. 

l<a,6<fc \ \ia + b - 1 



(B.2) 



det ((. " ,))= det ff^ + 'rV;/' 

l<a,b<k \ \l„ + b - 1 J J l<a.b<k \ \l„ + b - I J J 



which can be seen to hold by using Pascal's identity. Nevertheless, Pascal's identity 
in the quantum case, is slightly more complicated and reads: 



(B.3) 



n - 


h I 




n 




n 


J^ 


- 1 




.J + 1. 


J. 



"'^^See for instance 1441 Eq. (16.114) for a nice presentation. 
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and we thus see that in the quantum case the same kind of simphfication can not 
be shown to hold as simply as in the classical setting. 

To prove Eq. (|4.33|) nevertheless, first recall that in the space of symmetric 
polynomials the change of basis between the elementary symmetric polvnomialJ^ 
e,- and the Schur polynomials S\ is given by ( _38 Eq. (3.5)): 



(B.4) Sa = det(eA;-i+j), 

where A = (Ai > . . . > is a partition and A' its conjugate partition. This is actu- 
ally nothing else than Giambelli's identity, written for the symmetric polynomials 
and not just for the dimensions of the associated representations. To compute from 
this, note that if one considers x = {l,q, . . . , q"~^), then one has (IMI §1-3 Ex. 1. 
p.44): 

(B.5) 5a(x) = 

where n(A) is defined as: 

(B.6) n(A)^^(*-l)A, 

i>i 

and satisfies the following useful formulae: 



n 

A ' 




(B.7) 5]c(x)=n(A')-n(A) 

for the content (|S3, §1 Ex.3, pll ), and another one for the hook-lengths f |38|. 
Ex.2, p.ll ): 

(B.8) ^/i(a;) =n(A)+n(A') + |A|. 

xe\ 

Let us come back to our computation and particularize Eq. ljB.4|) to x 
(1, g, . . . , g"~"^), that gives: 



(B.9) 



— det q 



(A'-i + 3)(A'-i + 3-l) 



n 

A'; - i + j_ 



or, introducing + a — 1 = A'^, 



(B.IO) 



n(A) 



n , / {ia+b-i)i,ia+b-2) n 
, — det o 2 
AJ V [ia + b-1 

This is still not quite what we want. To proceed further note that according to 
Eq. HB.3|) one has: 



(B.ll) 



n + 1 
J + 1 



j(3 + l) 

= q ^ 



n 
J + 1. 



+ q" q 



j(j-i) 



Therefore, by multiple linear combinations of columns one can write: 



13 



Elementary symmetric polynomials are special cases of Schur polynomials of 1-column dia- 
grams, or Br = Sa^ in our notations. 



22 



YACINE DOLIVET AND MIGUEL TIERZ 



(ia+b-l)(ia+6-2) 



n 

ia + b 



(B.12) det 
Now, for convenience extract a factor q 
(B.13) 



= det I q 



(in+b-l)(i„+i>-2) 



n + b-l 
ia + b-1 



in each line to get: 



det 



{ia+b-l)(ia+b-2) 



b-1 

-5-1 



= q^^ 



det I q 



(ti-l)(2i„+b-2) 



-b 



To proceed further note the following property of the (/-binomial coefficients: 

'n + b-1' 
ia + b -I ' 



(B.14) 



q 





Vb 






Vb' 








ia - 


^b_ 



+ - 1) 



Therefore, once we write, 
(B.15) q 2 
it is easy to see that 
(B.16) det U^^^^ 



n + b 
ia + b_ 



n + b-l 
ia + b -I 



q 



b(ia+b) 



n + b 
ia + b 



q ^3 2 det 



Collecting everything we thus have 



det 



'n + b-l 
ia + b- 1 
From which we finally obtain: 



-+«(A) 



n + b-l 
ia + b-1 



(B.17) 



Appendix C. Notation for Chern-Simons quantities 



We give here some information about the Chern-Simons quantities that appear 
in the text, mainly in For more information, see ^ and references therein. To 

understand the origin of other quantities in , one has to take into account the 
constructions of Seifert homology spheres from surgery. Seifert homology spheres 
can be constructed by performing surgery on a link C in S'^ with n+1 components, 
consisting on n parallel and unlinked unknots together with a single unknot whose 
linking number with each of the other n unknots is one. The surgery data are Pj/qj 
for the unlinked unknots, j = 1, ■ ■ ■ , n, and on the final component, pj is coprime 
to qj for all j = 1, • • • , n, and the pj's are pairwise coprime. After doing surgery, 
one obtains the Seifert space M ^ X{^, - ■ ■ , |^). This is rational homology sphere 
whose first homology group Hi(M, Z) has order where 

n n 

qj 



(C.l) 



H 



pY^f, and P=Y{p,. 



Another topological invariant that will enter the computation is the signature of £, 
which turns out to be: 



(C.2) 



'^('C) = ^sign 



sign 
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For n = 1,2, Scifcrt homology spheres reduce to lens spaces, and one has that 
L{p,q) ~ X{q/p). Forn = 3, we obtain the Brieskorn homology spheres I](pi,p2,?33) 
(in this case the manifold is independent of 91, 92, 93)- In particular, S(2, 3,5) is 
the Poincare homology sphere. Finally, the Seifert manifold X{^, ^^^_^".^y^ , ^^5^), 
with m odd, can be obtained by integer surgery on a (2, to) torus knot with framing 
t. Note that in (|1.3|l the weight and root lattices of G are denoted by and Aj-, 
respectively. 

Finally, there is a phase factor in H1.3|l that comes from the framing correction, 
that guarantees that the resulting invariant is in the canonical framing for the 
three-manifold M . Its explicit expression is: 

(C.3) = 3sign(^-j +^12s(g„p,)-|, 

where cr(£) is again the signature of the linking matrix of C and s(p, q) is the 
Dedekind sum: 

(0.4) = i^g„.(=i)„.(I2?). 

n— 1 
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